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This Chapter gives a brief introduction to some basic aspects metals and superconductors in 
crystal without inversion symmetry. In a first part we analyze some normal state properties which 
arise through antisymmetric spin-orbit coupling existing in non-centrosymmetric materials and show 
its influence on the de Haas- van Alphen effect. For the superconducting phase we introduce a multi- 
band formulation which naturally arises due the spin splitting of the bands by spin-orbit coupling. 
It will then be shown how the states can be symmetry classified and their relation to the original 
classification in even-parity spin-singlet and odd-parity spin-triplet pairing states. The general 
Ginzburg-Landau functional will be derived and applied to the nucleation of superconductivity in a 
magnetic field. It will be shown that magneto-electric effects can modify the standard paramagnetic 
limiting behavior drastically. 

PACS numbers: 



I. INTRODUCTION 



Motivated by the discovery of the non-centrosymmetric heavy Fermion superconductor CePtsSii, the physics of 
unconventional superconductivity in materials without inversion symmetry has recently become a subject of growing 
interest. The lack of inversion symmetry, a key symmetry for Cooper pairing, combined with unconventional pairing 
symmetry is responsible for a number of intriguing novel properties. In a short time the list of new superconductors 
in this class has been enlarged by compounds such as Ulr— , CeRhSi^, CelrSi^, Y 2 C^ and Li 2 (Pdi_. r Pt 3 ;)3B&. In all 
listed heavy fermion compounds superconductivity appears in combination with a magnetic quantum phase transition 
suggesting the presence of strong electron correlation effects. Thus, it is widely believed that magnetic fluctuations 
are likely responsible for inducing here unconventional Cooper pairing. For other materials correlation effects seem 
to be less relevant. Nevertheless, some of them show unexpectedly features of unconventional pairing. While in some 
cases experimental results give rather clear suggestions on the gap symmetry, the definite identification of pairing 
states is far from concluded. 

The microscopic theory of superconductivity in metals without inversion has a long history predating these re- 
cent experimental developments 8 -^ ' 10 ' 11 . Specific aspects such as the possibilities of inhomogeneous superconducting 
state s 12 ' 13 and the magneto-electric effec t 14 ' 15 in this type of materials have been discussed already in the nineties of 
the last century. Moreover, general symmetry aspects of non-centrosymmetric superconductors have been addressed 
rather early in Refsi 16 ' 17 ' 18 . A wide variety of physical phenomena connected with non-centrosymmctricity have since 
been studied by many groups: 

• paramagnetic limitations of superconductivity and the helical vortex stat o 19 ' 20 ' 21 ' 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 ; 

• paramagnetic susceptibilit y 30 ' 31 ' 32 ' 62 and the magnetic field induced superconducting gap structure^ 3 -; 

• Josephson and quasiparticle tunneling— — . surface bound state s 36 ' 37 , and vortex bound states^ 8 -; 

• London penetration depth^ 9 . and the magnetic field distribution^ -; 

• effects of impuritie s 41 ' 42 ' 43 ; 

• upper critical fiel d 44 ' 45 ; 

• nuclear magnetic relaxation rat o 46 ' 47 ; 

• general forms of pairing interaction^ 8 -; 

• inhomogeneous superconducting states in the absence of external field 49 . 



In this Chapter we give a brief introduction to several topics in this context, leaving most of the special aspects of 
non-centrosymmetric superconductivity to other Chapters. 
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II. NORMAL STATE 

The absence of inversion symmetry is imprinted into the electronic structure through spin-orbit coupling effects. 
Already the normal state of non-centrosymmetric metals bears intriguing features which result from the specific form 
of spin-orbit coupling. In this section we discuss the electronic spectrum. 

A. Electronic states in non-centrosymmetric metals 

Our starting point is the following Hamiltonian of non-interacting electrons in a crystal without inversion center: 

Ho=Y, X [£( k )^+7(k)-o- a /3]4 a a k /3 (1) 

k a/3=U 

where at («ka) creates (annihilates) an electronic state |ka). Furthermore, £(k) = e(k) — /i denotes the spin- 
independent part of the spectrum measured relative to the chemical potential fi, a, f3 —] , J, arc spin indices and <x are 
the Pauli matrices. The sum over k is restricted to the first Brillouin zone. The second term in Eq. (fTJ) describes the 
antisymmetric spin-orbit (SO) coupling whose form depends on the specific non-centrosymmetric crystal structur o 50 ' 51 . 
The pseudovector 7(k) satisfies k) = — *y(k) and gj(g~ 1 'k) = 7(k), where g is any symmetry operation in the 
generating point group Q of the crystal (see below) . The usual symmetric spin-orbit coupling which is present also in 
centrosymmetric crystals yields a new spinor basis (pseudospinor) a, (3 in Eq. (fTJ), which retains the ordinary spin- 1/2 
structure with complete SU(2)-symmetry. This is different for the antisymmetric spin-orbit coupling. The effect of 
the antisymmetric spin-orbit coupling is a spin splitting of the band energy with k-dependent spin quantization axis 
which removes the SU(2)-symmetry. 

Depending on the purpose it is more convenient to express the Hamiltonian (JT|) in the initial 2x2 matrix form 
{spinor representation) or in its diagonal form (band representation) . The energy bands are given by 

£±(k)=£(k)±| 7 (k)| (2) 

with the Hamiltonian 

#o = X X ^(k)4 A c k A, (3) 

k A=± 

where the two sets of electronic operators are connected by a unitary transformation, 

&ka = X u aA(k)c kA , (4) 

A 

with 

l /, s /, \ \ (hl+A7^ Klx+ily)) 

(u T A(k), »u(k)) = , . ... i v — ■ (5) 

V2|7|(|7| +A 7z ) 

The normal-state electron Green's functions in the spinor representation can be written as 

G(k, u n ) = n A (k)G A (k, uj n ), (6) 

A=± 

where 

nx(k) = (7) 

are the band projection operators and 7 = 7/I7I. The Green's functions in the band representation have then the 
simple form 

G x (k,u n ) = -. l — — , (8) 

where w n — TtT(2n + 1) is the Matsubara frequency. 
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The Fermi surfaces denned by the equations £±(k) = are split, except at specific points or lines where |"y(k)| = 
is satisfied. The band dispersion functions £\(k) are invariant with respect to all operations of Q and the time reversal 
operations K = i^Ko (K$ is the complex conjugation). The states |k, A) and JsT|k, A) belonging to the band energies 
£x(k) and k), respectively, are degenerate, since the time reversal operation yields K\k, A) = (k) | — k, A), where 
*A(k) = — k) is a nontrivial phase facto r 11 ' 17 . For the eigenstates of H , defined by (O, this phase factor takes 
the form, 

= _ a7 (k)-» 7y (k) 

v / 7 2(k)+ 7 2(k) 

Finally we turn to the basic form of the antisymmetric spin-orbit coupling as it results from the non-centrosymmetric 
crystal structures. Here we ignore the Brillouin zone structure and use only the expansion for small momenta k leading 
to basis functions satisfying the basic symmetry requirements of 7 (k). For the cubic group Q = O, the point group 
of Li2(Pdi_2;,Pt x )3B, the simplest form compatible with symmetry requirements is 

7 (k) = 7o k, (10) 

where 7 o is a constant. For point groups containing improper elements, i.e. reflections and rotation-reflections, 
expressions become more complicated. The full tctrahcdral group Q = T^, which is relevant for Y2C3 and possibly 
KOS2O6, the expansion of 7 (k) starts with third order in the momentum, 

7 (k) = J [k x (k 2 y ~ k 2 z )x + k y {kl - k 2 x )y + k z {k 2 x - k 2 y )z]. (11) 

This is sometimes called Dresselhaus spin-orbit coupling 5 ^, and was originally discussed for bulk semiconductors of 
zinc-blend structure. 

The tetragonal point group Q = C^ v , relevant for CePtsSi, CeRhSi3 and CeIrSi3, yields the antisymmetric spin-orbit 
coupling 

7(k) = j±(k y x - k x y) + -f\\k x k y k z (kl - k%)z. (12) 

In the purely two-dimensional case, setting 7|| = one recovers the Rashba interaction^ which is often used to 
describe the effects of the absence of mirror symmetry in semiconductor quantum wells. 



B. de Haas - van Alphen effect 

An experimental way of observing the spin-splitting of the Fermi surface is the de Haas - van Alphen effect which 
can help to estimate the magnitude of the antisymmetric spin-orbit coupling^. The single-electron Hamiltonian ([1]) 
can be extended to include the magnetic field as follows: 

- ^Hir^la^ak/j. (13) 

The last term describes the Zeeman interaction for an external magnetic field H, with hb being the Bohr magneton. 
The orbital effect of the field can be included by replacing k — > k + (e/hc)A(r ),— where f = iVk is the position 
operator in the k-representation. 

The eigenvalues of the Hamiltonian (fTB")) are 

e A (k,H)=£(k) + A| 7 (k)-/i B H|. (14) 

There are two Fermi surfaces determined by the equations 

G(k,H)=0. (15) 

For certain directions and magnitudes of H there may be accidental degeneracies of the Fermi surfaces, determined 
by the equation 7(k) = /i^Ii. However, there are no symmetry reasons for such intersections. 

An important property of the Fermi surfaces (fTS"]) is the fact that their shapes depend on the magnetic field in 
a characteristic way, which can be directly probed by dHvA experiments. Note that while at H = time reversal 
symmetry guarantees (,\(— k) = Ov(k), the loss of time reversal symmetry for H ^ yields, in general, £a(— k, H) 7^ 
^A(k, H), i.e. the Fermi surfaces do not have inversion symmetry. 
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Including now the coupling of the magnetic field to the orbital motion of the electrons we derive in quasi-classical 
approximation the Lifshitz-Onsager quantization rules 5 ^ : 

Sx(e,k H ) = ^^[n + ax(T)}. (16) 
he 

Here S\ is the area of the quasi-classical orbit T, in the k-space defined by the intersection of the constant-energy 
surface e\(k) = e with the plane kh = kn (h = H/H). Moreover, n is an integer number (n> 1), and < chaIT) < 1 
is connected with the Berry phase of the electron as it moves along r 54 i 55 . The value of (^(r) does not affect the 
dHvA frequency discussed below. 

The dHvA signal contains contributions from both bands and can be approximately decomposed into the form, 

M osc = cos + A ) , (17) 

where A\ and (f>\ are the amplitudes and the phases of the oscillations. The amplitudes are given by the standard 
Lifshitz-Kosevich formula and the dHvA frequencies Fx are related to the extremal (with respect to kn) cross-section 
areas of the two Fermi surfaces, 

he 

Fx = —ST'- (18) 

In addition to the fundamental harmonics (|17p . the observed dHvA signal also contains higher harmonics with fre- 
quencies given by multiple integers of Fx- 

It is interesting to consider the field dependence of the band energies (fl4l , which yield 

Sl xt (U) = Sf*(0) + A x {h)H + B\(h)H 2 + ... . (19) 

Inserting this in Eq. (|17[) the term linear in H contributes to the phase shift, similar to the paramagnetic splitting of 
Fermi surfaces in centrosymmetric metals. The quadratic term is responsible of the magnetic field dependence of the 
dHvA frequencies. This is a specific feature of non-centrosymmetric metals which could be observable, if the Zeeman 
energy is at most of comparable magnitude as the spin-orbit coupling. 

fc2 ^,2 ^2 ^,2 

For illustration, let us look at the example of a three-dimensional elliptic Fermi surface with £(k) = 2m ^ + -^f- ~ e F, 
where m± and m z are the effective masses. The extremal (maximum) cross-sections of the Fermi surfaces (TT 
correspond to k z = 0. Introducing the Fermi wave vector fcp via cf — h 2 k F /2m±, we obtain 



ST*(H) = 7TA& 



(20) 



In this approximation we assumed that the Zeeman energy is small compared to the spin-orbit band splitting, which 
in turn is much smaller than the Fermi energy: [IbH <C |7_l|&f "C ep. Based on this result it is also possible to obtain 
an estimate of the strength of the spin-orbit coupling. 

We use the expressions (|18l) and ([20]) to calculate the difference of the dHvA frequencies for the split bands: 

For example, from the frequencies of the "a" and " j3" dHvA frequency branches in LaPtaSi^, F a = 1.10 x 10 8 Oe 
and Fj3 = 8.41 x 10 7 Oe, and m± ~ 1.5m, we obtain for the spin-orbit splitting of the Fermi surfaces: |7j_|fci? ~ 10 3 K, 
which is in reasonable agreement with the results of band structure calculation s 16 ! 56 . According to Eq. (|21|> . the 
magnetic field effect on F_ — F + in the range of fields used in Ref^ 6 - (up to 17T) should be of the order of a few 
percent. In this way the interplay of the Zeeman splitting and the spin-orbit coupling results in a deformation of 
the Fermi surface is responsible for a field dependence of the dHvA frequencies, an effect absent in centrosymmetric 
metals. 



III. SUPERCONDUCTING STATE 



In this section we turn to the discussion of some novel aspects of the superconducting state in non-centrosymmetric 
materials. Here we can consider only a few examples, while a wider range of other phenomena will be discussed in 
other Chapters of this book. 



5 



A. Basic equations 



After our introductory discussion of single-electron properties we now include electron-electron interactions to 
examine the implications of non-centrosymmetricity on Cooper pairing. Therefore we retain among all interactions 
only those terms corresponding to the Cooper channel and formulate it in the band representation. The general form 
is given by 

H int=2y^2 E T/ A 1 A 2 A 3 A 4 (k,k')cj c Ai C^ k A2 C_ k / !A3 Ck',A4- (22) 
k,k' Ai A2A3A4 

It is reasonable to assume Ai = A2 and A3 = A4, such that only intra-band pairing is considered. Inter-band pairing 
is usually suppressed, since the spin-orbit coupling induced band splitting would require that electrons far from the 
Fermi surfaces would have to pair which is unlikely, if the energy scale of the band splitting strongly exceeds the 
superconducting energy scale 57 . Introducing the notation A x = A 2 = A and A 3 = A4 = A', we obtain: 

H int = E Z)Vxv(k,k / )cJ e)A ci kiA c_ k ,,v<v 1 v, (23) 

kk'q XX' 

where 

V xx ,(k,k') = t x (k)tl,(k')V X y(k,k'). (24) 

Since under time reversal the creation and annihilation operators behave as 

Kc tx = *A(k)c! kiA , Kc k . x = i A *(k)c_ k ,A, (25) 

V\\'(k, k') represents the pairing interaction between time-reversed states. The amplitude V AA '(k, k') is even in 
both k and k' due to the anticommutation of fcrmionic operators and is invariant under the point group operations: 
Vxx 1 (5k, gk!) — V AA /(k, k'). The gap functions of the superconducting state can be expressed as, A A (k) = i A (k)A A (k), 
in each band, where A A transforms according to one of the irreducible representations of the crystal point grou p 59 ' 60 . 
The Gor'kov equations in each band read 

(iuj n -^ x (k))Gx(k,LJ n )+A kX Fl(k,LJ n ) = 1 (26) 
(uu n + U-k)) Fl(k,Lu n ) + Ai x G x (k,Lu n ) = 0. (27) 

The gap functions obey the self-consistency equations 

A k A = E E V\x (k, k') F\i (k', w „). (28) 

n k' A' 

The resulting Green's functions are then, 

G X (k,u n ) = ? "" +eA( ' k) ^-r- (29) 

(iu„ - 6 (k) ) {iw n + (-k)) - A kA A kA 

F A (k,w n ) = = ^r~. (30) 

(<«„ - UV)(iLu n + e A (-k)) - AkA A kA 

and the quasiparticle excitation energies for each band have the form 



Bu .^ »y tiiii , (31) 
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which becomes in case of time reversal symmetry, 

£ k A = ± V / a(k) 2 +A k AA kA . (32) 

These forms are analogous to that of a multi-band superconductor—, apart from the fact that in the non- 
centrosymmetric case the two bands to not possess spin degeneracy but rather correspond to a type of spinless 
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fermions, since their spinors on each band are subject to a momentum dependent projection. This distinction be- 
comes more apparent, if we write the Gor'kov equations in the initial spinor basis (spin up and down), 

(iw n - £(k) - 7k <x) G(k, «„) + A k #t(k, wn) = i (33) 

(^„+e(-k)+7-k^ t )^ t (k,^n) + Aj c G(k,w I1 ) = 0, (34) 

where £(k) = £ k — /i, 

G(k, u n ) = n+G+(k, c^„) + n_G_ (k, w B ), (35) 

Ft(k,w„) = g* {fl+Fl (k,w n ) + n_Fl(k,o; n )} , (36) 

A k = {n + A k , + + n_A k ,_}.g, (37) 

with 5 = i<7y. Examining the form of the gap function reveals that in the non-centrosymmetric superconductor both 
even-parity spin-singlet and odd-parity spin-triplet pairing are mixed, since no symmetry is available to distinguish 
the two. Therefore, we may write 

A A k , + + A k ,_ „ A k , + - A k ,_ „ . 

The odd-parity component is represented by a vector which is oriented along the vector / y k of the spin-orbit coupling. 
In this discussion the only approximation entering so far is the absence of inter-band pairing. 

B. Critical temperature 

For the discussion of the instability condition at the critical temperature and the topology of the quasiparticle gap 
we can use the symmetry properties of the pairing interaction matrix Vxv(k, k') mentioned earlier. The momentum 
dependence of the matrix elements can be represented in a spectral form decomposed in products of the basis functions 
of irreducible representations of Q. It is generally sufficient to consider only the part of the pairing potential based 
on one irreducible representation T corresponding to the superconducting state with maximal critical temperature^. 
In a simplified formulation this could be represented by even basis functions on the two bands, 0+,j(k) and 0_.j(k), 

d 

V xx ,(k,k') = -^^^(k^^k'), (39) 

i=l 

While (f>+,i (k) and </>_.j(k) both belong to same symmetry representation, their momentum dependence does not 
have to be exactly the same. The basis functions are assumed to satisfy the following orthogonality conditions: 
(</>^ j(k)^A,j '(k)) a = where the angular brackets denote the averaging over the A-th Fermi surface. The coupling 
constants V\\> form a Hermitian matrix, which becomes real symmetric, if the basis functions are real. The gap 
functions take the form 

A A (k) = ^r ?A)i A)i (k), (40) 

i=l 

and -q\ t i are the superconducting order parameter components in the A-th band. 

As an example, consider a superconducting state with the order parameter transforming according to a one- 
dimensional representation A A (k) = i]x4>\(^-)- The linearized gap equations Eq. (|28|) acquire simple algebraic form 

V+ = (9++V+ + 9+-V-)Si{T), 

V- = (g-+V++9—V-)Si(T), (41) 

where 



.9A A' = V\\> Nq\i , 



(42) 
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and No\ — (\4>\(]t)\ 2 No\(k)) \ is the weighted average angular dependent density of states over the A-th Fermi surface. 
Note that for multidimensional representations (dimensional dp) due to the crystal point symmetry the values of 

Nox = (|0A,i(k)| 2 A r O A(k))A are equal for all components i = 1, ..dr and all components 77^,1, ■ ■■r]\,dr separately satisfy 
the same system of equations (|4Tj) . 
The function Si(T) is 

Sl (T) = 2nTj2-=l^, (43) 

where In 7 = 0, 577 ... is the Euler constant, e c is an energy cutoff for the pairing interaction, which we assume to be 
the same for both bands. From Eq. (|41|) we obtain then the following expression for the critical temperature: 

T c0 = ^e Xp f-iY (44) 
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where 




g= o ~ +\ * n " +9+-9- + (45) 



is the effective coupling constant. For multidimensional representations the critical temperature is the same for all 
dr components of T)\j of the order parameter. The particular combination of amplitudes r\\.i in the superconducting 
state below T c is determined by the nonlinear terms in the free energy or self-consistent equation, which depend on 
the symmetry of the dominant pairing channel. 

The solution of Eq. (|41[) (r]+, r]J) corresponding to the eigenvalue S\{T C ) determines two unequal order parameter 
components A^(k) = ^^(k). In the spinor representation (|3"5|) both singlet and triplet parts of the order parameter 
are present. Pure singlet or pure triplet pairing occurs only under rather restrictive conditions. First, the momentum 

dependence of the gap function in both bands is the same Ax(k) = r)\(j)(k.). Second, g++ = g and g^ = g |_ is 

realized. Then we obtain two solutions of equations Eq. (|4"Tj) with 



V+ = V- . (46) 



V+ = -V- ( 47 ) 

The critical temperature of the state (|4€>[) corresponding to the singlet part is T c s = (2e c /7r)e _1 / 9s , where g s — 

g ++ + <7_| The critical temperature of the spin triplet state (|T7]) is T c ' = (2e c /n)e^ 1 ^ 9t , where gt = g++ — g^ If 

<JL| > then T c s > T* and the phase transition occurs to the state (|46|) . While at g^ < we see that T c ' > T c s 

and the phase transition occurs to the state (|4"7| . 



C. Zeros in the quasiparticle gap 



On the one hand, the zeros in the gap for elementary excitations are dictated by the symmetry of superconducting 
state or its superconducting class which is a subgroup TL of the group of symmetry of the normal state Q x K, x f/(l), 
where Q is the point group, K, is the group of time reversal, U(l) is the gauge group. The procedure to find 
symmetry dictated nodes is described in Ref.— . Let us consider the possible superconducting states (|40|) and their 
nodes for CePt^Si with point group symmetry C^ v . This group has four one-dimensional irreducible representations, 
Ai, A2, Bi, B2, and one two-dimensional representation, E. Examples of even basis functions of these irreducible 
representations are 



r 


<Mk) 


nodes 


At 


e x + e y + ck 2 z 




A 2 


k x k y (k 2 — ky) 


k x — 0, ky = 0, k x = i/^y 




h 2 h 2 

™x y 


kx — -^-ky 


B 2 


kx ky 


k x — 0, ky = 
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For E state the basis functions are </>Ei(k) = k x k z and 0£2(k) — k y k z and it leads to the order parameter for the 
Fermi surfaces A, 

A A (k) = ti\i<P\ei (k) + 77A,20A,s 2 (k). (48) 

The symmetry of superconducting state and the corresponding node positions depend on the particular choice of 
amplitudes r)\.i. Under weak coupling conditions the combination generating least nodes is most stable, corresponding 
here to (t]\,i, T]x,2) = ^a(1,±*) and is a time reversal symmetry violating phase with a line node on the plane k z = 
and point nodes at k x = k y = 0. 

On the other hand, the fact that their gaps on the two Fermi surfaces are composed of an even and an odd parity 
part, can also lead to nodes which are not symmetry protected, as discussed in Ref.— . 



D. The amplitude of singlet and triplet pairing states 

The coupling constants V X y we have used in previous considerations can be expressed through the real physical 
interactions between the electrons naturally introduced in the initial spinor basis where BCS type Hamiltonian has 
the following form 59 

kk'qa/3 7 5 

+y^(k,kO(jCT l CT 2 ) Q/3 (icr J CT 2 )t (5 ]aj i+qQ a 1 l k)g a_ k ^ 7 a k ' +q ^, 

here the amplitudes V 9 (k, k') and V^"(k, k') are even and odd with respect to their arguments correspondingly. The 
unitary transformation (01 transforms the pairing Hamiltonian (|49[) to the band representation (f2"2"| . If we neglect 
inter-band pairing, it is reduced to ()23|1 and (|24|) with the amplitudes given by the following expression 



1/ A v(k,k') = ^ S (k,k')(<7 + a x ) X x> + ^(k,k')7i(k)7j(k')(<7 - a x )xx>- (50) 



The explicit derivation is given ini&, where the similar procedure was also made for more general interactions mediated 
by phonons or spin fluctuations. It can be shown that the pairing given by the amplitude V 9 (k, k') in the initial 
spinor basis including the simple s-wave pairing V 9 (h, k') = const does not induce any inter-band pairing channel. 

To illustrate the origination of the singlet and triplet pairing channels let us consider a superconductor with 
tetragonal symmetry C^ v and Rashba spin-orbital coupling 7(k) = 7i(z x k), for a spherical Fermi surface. We 
describe the pairing by the following model, which is compatible with all symmetry requirements: 

V(&k!) = -V a , 

^(k,k') = -K(7 i (k)7i(k / )), (51) 

where V g and V u are constants. This type of pairing interaction yields the superconducting state with full symmetry 
of the tetragonal group C^ v transforming according to unit representation A\ both in singlet and in triplet channels. 
With Eq. ([ST)]) we arrive at the band representation: 



V3, A '(k,k') = -lv g {a + <t x ) xx > - ~V u (a -a x )xx', (52) 



Thus, this pairing interaction is even simpler than that (eqn. (|39l) ) considered in the previous subsection. So, in our 
model the gap functions in the two bands Eqs. (|40[) are: A^(k) = rjxLp^ (k) wth k independent functions ipA t (k) = 1. 
The amplitudes rj\ satisfy the equations 

Vx = Y] 9xx r ttT V —f= ? X \ , (53) 

y n V^l + Vx> 



where 

g±± = g± T = 1J h r ^ N 0T- (54) 



Vg + V u ^ _V g - V u 



The critical temperature is given by Eq. (|44[) . 
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According to the Eq. (|38|) the singlet and triplet parts of the order parameter are determined by the order parameter 
amplitudes in different bands. For the ratio of triplet to singlet amplitude in vicinity of T c we find: 

r - - 7 ?- _ 2 -9+- + £±± - a— - VP (55) 

V+ + V- 2.g + _+.g__-.g ++ + VP' 

where T> = (g++ — g ) 2 + 4g_| g It is easy to see that for V u = the triplet component of the order parameter 

vanishes identically: r = 0. On the other hand, for V g = the singlet component of the order parameter disappears: 
= 0. Generally the relative weight of singlet and triplet component in the order parameter depends on the ratio 
of pairing interactions decomposed into even and odd parity channel. 
A simple BCS type of model with 

V a (k,k') = -V g and V£(k,k') = (56) 

yields in the band representation 

V x B x ? s (k,k>)=-±V g (a +a x ) xx , (57) 
and gives rise to purely spin singlet pairing within our notion. 



E. Ginzburg-Landau formulation 



The Ginzburg-Landau theory is a very efficient tool to discuss a wide variety of phenomena of the superconducting 
state, in particular, the instability conditions at the critical temperature. We will derive the Ginzburg-Landau 
functional from a microscopic starting point with the aim to address in the following chapter the influence of the 
magneto-electric effect on the nucleation of superconductivity in a magnetic field, i.e. the modification of paramagnetic 
limiting in a non-centrosymmetric metal. 

For this purpose we extend the self-consistent equation (|28|) to the case where magnetic fields are present and the 
superconducting order parameter has a weak spatial dependence, 

A A (k, q) = T]T £ V» ( k ' k ') G "( k '- G 4-k' + q, -U} n )A v (k', q ). ( 58 ) 

n k' v 

Near the critical temperature where one can use the normal metal Green functions G x (k,u) n ) which yield then the 
linearized gap equation to examine the instability condition. This equation can be derived from the free energy 
functional of the form 



(59) 

d 3 k 



/d k 
^ 5 A*(k,q)(? v (k J w„)G v (-k + q l -w n )A I/ (k,q) 



The corresponding normal metal electron Green function G\(k.,u) n ) = (iuj — £\(k, H)) 1 in a magnetic field is deter- 
mined by the electron energies (fT4l) . 



£ A (k,H) = £(k) + A| 7 (k) - mbH| « ^(k) - m A (k)H, (60) 

where ^ A (k) = £(k) + A|'y(k)| and the second term on the right-hand side is the analog of the Zceman interaction for 
non-degenerate bands^ with the form: 

m A (k) = A MB 7(k), (61) 

which is valid everywhere except for the vicinity of the band crossing points, where the approximation of independent 
non-degenerate bands fails. In the standard centrosymmetric metal the magnetic field splits the Fermi surfaces into 
majority and minority spin surfaces. Here the spin-splitting is imposed at the outset by the spin-orbit coupling. The 
effect of the magnetic field is a non-centric deformation of the band and shape of the Fermi surfaces. As we will see 
this will influence the superconducting condensate nucleated in a magnetic field. 
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The normal electron Green function is then approximated as 

Gx(k,w„) = , - 7 r^- j—. (62) 

iuj n - U(k) + mA(k)H 

Since the gap function weakly depends on energy in the vicinity of the Fermi surface, one can integrate the products 
of two Green's functions with respect to £\ = £\(k): 

N \ J d£_ x Gx(k,uj n )Gx{-'k + q, -ui n ) = nN 0X Lx(k, q,w„), (63) 



where 



depends only on k, the direction of k, 



L x (k,q,uj n ) = j— — — — (64) 

\u n \ +iilx[k, q) signw„ 



r> A (k,q) = ^&-m A (k)H, (65) 

with v A (k) = d£x{k)/dk being the Fermi velocity in the A-th band. 

The Ginzburg-Landau free energy in usual coordinate representation (as well as the Ginzburg-Landau equations) 
can be obtained from the Taylor expansion of Eqs. and (|BT)|) in powers of fi A (k, q), by the replacement 

q^D = -iV r + 2eA(r) (66) 

in the final expressions. In the following we will put h = c = 1. The special form of f2 A (k, q) introduces a novel 
gradient terms in the free energy of non-centrosymmetric superconductors. Instead of powers of q it contains powers 
of f2 A (k, q). This can can lead to the formation of nonuniform superconducting state known as helical phases and 
magneto- electric effect in a magnetic field. 



IV. MAGNETO-ELECTRIC EFFECT AND THE UPPER CRITICAL FIELD 



The term "magneto-electric effect" in non-centrosymmetric superconductors encompasses several intriguing features. 
It has been discussed on a phcnomcnological level by introducing additional linear gradients terms to the Ginzburg- 
Landau free energy, so-called Lifshitz invariants, like 

^(^KijHiD^ir) (67) 

Here ?7(r) denotes the order parameter of superconductor, H is magnetic field and D = — iV — 2eA is the gauge- 
invariant gradient. First predicted by Levitov, Nazarov and Eliashberg 9 -, the magneto-electric effect was studied 
microscopically by several author o 10 ' 14 i 15 ' 19 i 23 . In this context several observable effects have been predicted: (i) the 
existence of a helically twisted superconducting order parameter in a magnetic field in two and three dimensional 
cases and spontaneous supercurrents in a 2D geometr y 15 i 19 i 21 i 23 i 26 i 29 and near the superconductor surface^ as well 
as along junctions of two superconductors with opposite directions of polarization^!, (ii) the augmentation of the 
upper critical field oriented perpendicular to the direction of the space parity breaking 23 ' 2 ^, (iii) magnetic interference 
patterns of the Josephson critical current for a magnetic field applied perpendicular to the junction^ 6 .. 

The presence of Lifshitz invariants (|67[) . however, can mislead to invalid conclusions and a careful analysis of 
different contributions to an effect is mandatory. This is indeed true for the influence of the magneto-electric effect 
on paramagnetic limiting. Moreover, the notion of helical phase has to be considered carefully as it may imply wrong 
pictures. In this section we would like to give insight into this subtleties by discussing the magnetic field dependence 
of the effective critical temperature in the Ginzburg-Landau framework. 



A. One band case 



Before considering the intrinsic multi-band situation due to the spin splitting of the electron band, for simplicity 
we restrict ourselves to a one-band situation, i.e. we ignore one of the two bands. This band shall be characterized 
by an isotropic density of states at the Fermi energy, iVo+(k) — N+. 
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The Ginzburg-Landau free energy for this one-band case with a one-component order parameter can be derived 
from Eqn. 1(55]) 

F=\J ^|^--7Vo + 5 1 (T)+iVo + 53(r)<(0 2 (k)f7(k,q)) 2 )|| 7? (q)| 2 , (68) 

where we restrict to the second order term a 23 i 43 . This is sufficient to analyze the instability conditions. Here, <fi(k) 
describes the superconducting state and is an even function belonging to one of one-dimensional representations of 
the point group of the crystal, f2 is given by (|65j) . (...) means the averaging over the Fermi surface, the function Si(T) 
is given by eqn. (|43|) and 

n ' 1 

The Ginzburg-Landau free energy functional in real space can be obtained through a Fourier transformation and leads 
to 

F = Jd 3 r {a{T - T c0 )\ V \ 2 + V * [K X (D 2 X + D 2 y ) + K 2 D 2 Z + K^H.Dj + (1,11,11, r,} , (70) 
where a = N Q+ /2T c0 , T cQ = (2 7 e /7r) exp (-2/y++iV 0+ ), 

Ki = ^(0 2 (k), 2 (k)), K 2 = ^(^(k)^(k)), (71) 

Kii = - » BN » Sa {?m<Mvi<n), Qn = &^(4?<mmto)- w 

The term linear in H incorporates the magneto-electric effects while the term quadratic in H describes the paramag- 
netic effect. Therefore Qy|?j| is connected with the change of the paramagnetic susceptibility in the superconducting 
phase compared with the normal state (Pauli) susceptibility. In particular, Qij vanishes when there is no change 
of the paramagnetic susceptibility. These coefficients have to be compared with those of a spin singlet state of a 
centrosymmetric superconductor, = (j,gNoSs/2 which, assuming Nq + = No, is larger than above Qij due to the 
fact that 1 = (<^ 2 (k)) > (</> 2 (k) 7i (k) 7i (k)). 

We consider now the two illustrative cases, the point group C^ v and D4 which are characterized by the pseudovectors 

7(k) = J±(k y x - k x y) + ^\\k x k y k z (k 2 x - k 2 )z for C iv , 
7(k) = j±(k x x + k v y)+j\\k z z for D 4 . 

For symmetry arguments and using above expressions we find the following relations for the coefficients, 

K xy = —K yx 7^ and Kij = otherwise 

Qxx = Q yy 7^ Qzz > and Qij = otherwise 

for C± v where presumably \K ZX \ <C \K xy \ and Q zz -C Q xx due to large number of nodes in the k-dependence of the 
7||-part of 7(k) , and 

K xx = K yy 7^ 0, K zz 7^ and = otherwise 



(75) 

Qxx — Qyy 7^ Qzz > and Qij = otherwise 



for D 4 



1. Symmetry Civ, H || z 

In the case of C± v for the field directed parallel to z-axis H = H(0, 0, 1) the terms linear in gradients and H is 
absent. The standard solution r\ = e lqyV f(x) of the GL equation 



cO) 



cP_ 
dx 2 



i — 

' dy 



2eHx 



Q ZZ H 2 \n = 0, 



(76) 
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is degenerate in respect to q y . The magnetic field dependence of critical temperature is 

Tc = T c0 - - Q^H 2 (77) 

a a 

Both the orbital (linear in H) and paramagnetic (quadratic in H) depairing effect are present. Compared to the 
ordinary spin-singlet case, however, the effect of the paramagnetic limiting is weaker here due to Q zz < Q z °). It is 
important to note here that no magneto-electric effect comes into play here. 



2. Symmetry D4, H || z 

The situation is quite different for the uniaxial crystals with point symmetry group D 4 (or D§). The GL equation 
includes gradient terms in the field direction and acquires the form 



a{T - T c0 ) + K 1 



dx 2 



d v 

-i— + 2eHx 
dy J 



r, rr 9 „ 8 2 

+iK zz H—-K 2 — 
oz oz l 



QzzH 2 \n = o. 



(78) 



The solution can be written as 



T] = e iq y y e iq * z f(x), 



(79) 



which remains degenerate with respect to the wavevector q y , but not with respect to q z which is used to maximize 
the critical temperature to 



T r — T C Q — 



2eK 1 



H 



a 



K 2 

ZZ 

4K 2 



H 2 



This corresponds to the finite wavevector 



Qz = 



K ZZ H 



2K 2 



(80) 



(81) 



Note that this wave vector could also be absorbed into the vector potential without changing the physically relevant 
results: A — > A + V% with \ = —Qzz/2e. 

The simple paramagnetic depairing effect is weakened due to magneto-electric response of the system. Adjusting the 
nucleation of the superconducting phase to the shifted Fermi surface, as incorporated in the wavevector q z , recovers 
some of the strength of the nucleating condensate. This is a specific effect of the non-centrosymmetric superconduc- 
tor and has its conceptional analogue in the FFLO phase for ccntrosymmctric spin singlet superconductors, where 
the condensate also nucleates with finite momentum Cooper pairs in order to optimize the pairing of degenerate 
quasiparticles on the split Fermi surface. 



3. Symmetry Cav, Hli 

Now we turn the magnetic field into the basal plane H = H(cos ip, sintp, 0), and impose a gauge to have the vector 
potential A = Hz(smf 7 — cos ip, 0). The corresponding GL equation take the form 



ja(T - T c0 ) + K x (p\ + D 2 y ) - K 2 ^ + K xy (H x D y - H y D x ) + Q XX H 2 ^ V = 0, 



where 



(82) 



d d 
D x = -%— + 2eH y z, D y = -%— - 2eH x z. (83) 

Like in ordinary superconductors the solution of this equation have the Abrikosov form 

r)(r) = exp[i(pxr)J/(z), (84) 



13 



where we write p = pH/H as a vector parallel to the magnetic field ((p x r) 2 denoting the z-component of the vector 
p x r), and f(z) satisfies the resulting renormalizcd harmonic oscillator equation 

I o(T - T c0 ) + K 1 (2eH) 2 {z - z Q ) 2 - K 2 -^ + t Q xx - ^ \ H 2 \ f(z) = 0, (85) 



with the shifted equilibrium position 



- = '2.//, 1 iP+^ H )- 



Thus, the vector p is absorbed into the shift zq and does not appear anywhere else in the equation. Then the 
corresponding eigenvalue determines the magnetic field dependence of optimized critical temperature: 

a \ 4a i / a 

In the used gauge the eigenstates are degenerate with respect to p and acquire the same structure as the usual Landau 
degeneracy. Nevertheless, the characteristics of the non-centrosymmetricity incorporated in the A^-terms appears in 
the expression of T c . Similar to the previous case of D4 with H || z the magneto-electric effect yields a reduction of 
the paramagnetic limiting term. This renormalization is surprisingly strong in general, as we can see when we return 
to the expressions which we had derived for the different coefficients. We obtain for the last term in Eq. (|87|) , 



~K 2 

xy 




H 2 


4Ki 


Qxx 


a 



<^(k)t£(k)> 



fi 2 B H 2 N+S 3 
2a 



Considering the simplified picture of a parabolic band with v(k) = k/m* and a Rashba spin-orbit coupling 7 = k x z. 
(setting 7 2 (k) = 0) we find the amazing result that the two terms exactly cancel and the paramagnetic effect is 
completely suppressed. This effect can be immediately obtained, if we perform the gauge transformation 

q ^ q+ 2 M Bm*(zxH) 

kp 

already in Eq. (|65p and so eliminating the paramagnetic term at the outset. However, it is important to notice that 
this exact cancellation is a consequence of the simplified forms of the band structure and the spin-orbit coupling 
term. Taking more realistic band structure effects into account it is obvious that this identity does not hold anymore 
in general. Nevertheless, our results suggests that the magneto-electric effect can, in principle, yield a substantial 
contribution to eliminate the paramagnetic limiting also for fields in the basal plane. 

4- Symmetry D4, H _L z 

It is easy to see that this case is analogue to the situation for the field along the z-axis and has only quantitative 
differences. Thus also here we encounter a reduction of the paramagnetic limit due to the magneto-electric effect 
yielding 

r ._ T- _2s£^ jr+ ('£i_ a J5!, (90) 



where also the same considerations concerning the gauge freedom apply as in the case of H || z apply. 

5. Two-dimensional case, symmetry Ci v , H _L z 

The simplest way to pass from 3D to 2D situation it is to introduce S(z) function potential well into 3D GL equation 
(|82|) . It is equivalent to the theory used by Tinkham^l for the calculation of the upper critical field in a thin film with 
thickness d << £ for a field parallel to the film. Thus, we consider the instability equation 

[a{T ~ T cQ ) - K ± (D 2 X + D 2 y ) - K 2 -§^ + K xy (H x D y - H y D x ) + Q XX H 2 

~^S(z)} V = Q, (91) 
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where d is a length of the order of the film thickness that is in pure 2D case it is an atomic scale length. This eigenvalue 
equation has the solution 



rj(r) = Aexp [i(p x r) z ] exp ( - — 



(92) 



where p = pH/H is a vector with arbitrary length directed along magnetic field. This then determines the critical 
temperature as a function of the applied magnetic field. 



a(T - f c0 ) + K x {2eHf{{z - z Q ) 2 ) + Q xx 



K 2 

x y 

1^7 



H = 0. 



(93) 



Here T c q is the critical temperature in the absence of a magnetic field, corresponding to d 2 — K2/a(T C Q — T c q). 
Moreover, (...) denotes the expectation value using the wave function exp(—\z\/d) and zo is determined by the same 
expression as in the 3D case 



zo = ^eH)- 1 
Hence, we obtain for the critical temperature 



2Ki 



H 



K 2 

xy 

4^7 



a 



^(z 2 + d 2 /2)(2eH) 2 . 



The critical temperature reaches obviously a maximal value at zo = 0, i.e. for 



P 



K 

2~K, 



H. 



(94) 



(95) 



(96) 



The upper critical field shows also here the square root temperature dependence usual for thin films in a parallel 
magnetic field 61 . Under special conditions (e.g. rotation symmetry around the normal vector of the film) the expression 
in the square parenthesis in Eq. (|95[) may vanish, as described above. Then, unlike in usual superconductors, non- 
centrosymmetric superconductors follow the standard Tinkham behavior unchanged by paramagnetic contributions. 

In view of strong inequality d <C l/V2eH the complete suppression of 2D superconducting state (T C (H) = 0) is 
reached in the field which exceeds the orbital critical field in the 3D case (1871). 



B. Two band case 



While the one-band picture discussed so far gives useful insights into the influence of the magneto-electric effect 
on the upper critical field, in particular, in the context of paramagnetic limiting, in reality there are at least two 
split bands whose Fermi surface allows for the nucleation of a condensate in a finite magnetic field. In the two-band 
picture the situation is somewhat more complex, so that we restrict here to a few aspects only which, we believe, 
are relevant in this context without attempting to give a complete overview. We base our analysis on the formalism 
introduced for the homogeneous superconducting phase in section 3.2. We use also an order parameter belonging 
to a one-dimensional representation on the two Fermi surfaces, labeled by A = ±, A>(k, r) = r)\(r)(f>\(k). Moreover 
we restrict our discussion of the case of the point group C^ v with an in-plane magnetic field. Then the linearized 
Ginzburg-Landau equation is given by 



n+ = 9++[Si{T) - L+]r] + +.g + _[5i(T) - L_]r?_, 
V- = 9-+[Si(T) - L + ] V+ +g__[S 1 (T) - L_]ry_, 



(97) 



with the operators L\, 



Lx = N, 



OA 



K lx {Dl + D 2 ) - K 2X 



^K X y\{H X Dy — HyD X ) + Q X x\H 2 



(98) 



where the coefficients are defined through the straightforward generalization of Eqs. ([7T|) and (f72|) to the two-band 
case with the gap functions <^(k), the Fermi velocity components Ux,i(k) and the densities of states N\ taken in the 
corresponding band. 
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Similar to the one-band case the solutions of this equation system can be cast into the Abrikosov form 

exp [i(p x r) z ] 



r?_(r) 



/+(*) 
/-(*) 



(99) 



where again p = pH/H and the functions f+(z), f-(z) satisfy to the system of equations 

/+ = g++[Si(T) - M+}f+ + g+- [Si(T] - M_]/_, 
/_ = g-+[Sx(T) - M+]f+ + P __[Sx(T) - M_]/_. 

Using the same gauge as in the one-band example, the new operator M\ is then 



M A = N, 



OA 



KiA(2 ej ff) 2 (z-z A o) 2 



^QxxX 



K , 



4 if. 



H 2 



z 0X = (2eH)- 1 (p+X^H 



(100) 

(101) 
(102) 



As in the one-band case the eigen states of this system possesses the Landau degeneracy represented through the 
equilibrium positions of the coupled harmonic oscillators, zq + and zo_, which both depend on p. Through the 
substitution 



z = Z 



P 
2eH 



we can formulate the equation system so that p is eliminated and zq\ 



z. 



OA- 



z 



OA 



X- 



(103) 



The general solution of Eq. (|100[) can be found only numerically. Here we limit ourselves to a variational solution 
of the form, 



f + {Z) = C + e*v{- 



f-(Z) = C_exp 



eH^/K^jZ - Z _ 



(104) 



In the following calculations, taking into account that the band splitting is much less than the Fermi energy |7_l|&f *C 
£p, we neglect the difference of the Fermi velocities v\(k) and the densities of states Nox of the two bands. This case 
the values of Zq\ = XZq for different bands differ each other only by sign. Returning to the free energy functional 
and integrating over Z we obtain new variational equations for the coefficients C+ and CL : 



C+ = g++[Si(T) - M(H)]C+ + Ig+- [Si(T) - M(H)]C„, 
C_ = Ig-+[Si(T) - M(H)]C++g^[S 1 (T) - M(H)]C. 



with 



and 



M(H) = JV" 1 



2eHy/K 1 K 2 + Q 2 



I = exp (-2eHZ$) 



K 2 

xy 



H 1 



(105) 



(106) 



(107) 



The system (|105j) has the same form as the system (|41|) in the absence of magnetic field. Hence for the critical 
temperature we obtain 



T c = f c0 (l-M) 



(108) 



where he temperature T c0 is given by the same formula (|44[) as T c0 taking in mind the substitutions g^ — » g^ = 

Ig^ , g |_ — > g |_ = Ig |_. The product 2eHZ 2 w eHra* 2 / (kprn) 2 ) is much less than unity generally, except for 

heavy fermion or layered superconductors. 
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Thus, in the two-band situation the paramagnetic suppression of superconducting state oc —Q XX H 2 is substantially 
weakened by the magneto-electric effect oc K^. H 2 /&K\. The latter has a finite value so long we work in the limit 
HbH < \j±\k F . 

This conclusion is qualitatively valid in general, although our discussion was done by a simple variational approach 

only. Moreover, assuming g ++ = g and = g |_, as it was done in the absence of magnetic field (see eqns. 

(|46|) and (|47|) ). we come to the solution of (|105|) with either pure singlet or with pure triplet pairing. The effect 
of paramagnetic limiting takes place identically in both cases. This underlines directly that the weakening of the 
paramagnetic limiting in the non-centrosymmetric superconductors is not connected to the formation of a mixed 
singlet-triplet state. The important point is the spin-orbital splitting of the bands. The Pauli spin susceptibility 
of the quasiparticles in whole space between the two Fermi surfaces is not changed in the superconducting state in 
comparison of its normal value that leads to the weakening of the paramagnetic suppression of the superconducting 
state. 

Various approximations to the two-band model have been used in literature, some of which can obscure the subtleties 
of non-centrosymmetric superconductors. In particular, the discussion of the spin susceptibility in the superconducting 
phase as discussed i n 30 i 31 ' 32 have to be considered with caution in view of the magneto-electric effects which are 
neglected there. The adjustment of the superconducting state to the field-induced shifts of the Fermi surface yield a 
correction to the spin susceptibility which is not negligible as our discussion in the Ginzburg-Landau regime show. 
The subtle two-band effects, however, often make quantitative predictions difficult—. 



V. CONCLUSION 



In this chapter we have given an overview on some theoretical aspects of non-centrosymmetric superconductors. 
Unlike symmetric spin-orbit coupling found in centrosymmetric metals, the antisymmetric spin-orbit coupling has 
a spectacular influence on the electronic bands through a specific spin splitting of the quasi-particle states. Su- 
perconductivity as a Fermi liquid instability is naturally influenced by such a modification of the electronic states. 
Superconductivity then has always a multi-band character. Moreover, parity does no longer provide a good quantum 
number to classify the superconducting phases. 

One of the physically most remarkable aspects of non-centrosymmetric superconductivity is connected with 
magneto-electricity, the peculiar connection between supercurrents and spin polarization. We have considered one 
aspect in this context, namely its influence on paramagnetic limiting. This effect is of interest in strongly correlated 
electron systems where the coherence length is generally small due to the enhanced masses like in heavy Fermion 
compounds. Here ordinary orbital depairing in a magnetic field is weak, such that the upper critical field reaches 
magnitudes where paramagnetic limiting through spin polarization becomes visible. Interestingly already on the basis 
of symmetry considerations it is possible to arrive at very interesting predictions which are borne out in some of the 
non-centrosymmetric heavy Fermion superconductors. 
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